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Abstract. In ground based infrared imaging a well-known technique to reduce the influence of thermal and 
background noise is chopping and nodding, where four different signals of the same object are recorded from 
which the object is reconstructed numerically. 

Since noise in the data can severely affect the reconstruction, regularization algorithms have to be implemented. 
In this paper we propose to combine iterative reconstruction algorithms with robust statistical methods. 
Moreover, we study the use of multiple chopped data sets with different chopping amplitudes and the according 
numerical reconstruction algorithm. Numerical simulations show robustness of the proposed methods with respect 
to noisy data. 
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1. Introduction 



Data collected in infrared ground based astronomy are 
affected by atmospheric and telescopic thermal back- 
ground noise. A common approach for noise reduction is 
chopping and nodding l)Emersonl Il993 iRobberto et~al] 
120051) . 



In cho pping the secondary mirror of the telescope is 
moved fcf.l Bertero et al.l ()2003a|) ) and signals are recorded 
for two different tilt angles. We denote the positions before 
and after tilting by A and B = A+h, and the recorded sig- 
nals by Si and 5*2, respectively. Afterwards, in a nodding 
procedure the telescope is moved to the position C = A—h 
and the signal S3 is recorded. Then a second chopping pro- 
cedure is applied and a signal 5'i is recorded at position 
A. 

In literature the vector h is referred to as the chop- 
ping amplitude. The data visualized after chopping and 
nodding is 

f — Si — S2 — S3 + Si 

and the image u to be reconstructed satisfies 
2u{.)-u{.-h)-u{. + h)^f{.) . 

When observing a point-like and isolated object the chop- 
ping amplitude can be chosen in such way that at positions 
B and C only empty sky is recorded in which case the nu- 
merical reconstruction algorithm is dispensable. However, 

to: Frank 



in practice the objects observed are often extended or 
other objects are nearby. In these cases the reconstruc- 
tion step is mandatory. 

Reconstructio n from chopped and nodded data was 
first discussed in iBecke rs (1994(), where a Fourier-based 
reconstruction has been suggested. 



In iBertero et all l)l999l I2OO0I l2003albl) iterative re- 
co nstruction a l gorith ms have been investigated. Recently, 
m et al.l l(200.ll) reconstruction algorithms based on 

wavelet decomposition have been proposed. 

Experiments have shown that noise can sev erely af- 
fect the reconstruction process fsee lKa,eufllll"995|) . In this 
paper we investigate two different approaches for robust 
reconstruction in the presence of a significant amount of 
noise. 

The outline of this paper is as follows: in Sect. [21 we 
derive the model of chopping and noddi ng and review the 
reconstr uction algorithms introduced bv lBertero et al.l (cf. 
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In this paper we use iterative regularization techniques 
as well and thus Bertero's work is the most relevant to 
compare with. However for the sake of noise robustness 
we combine iterative regularization techniques with robust 
filtering methods from statistics (cf. Sect.|31l. Moreover, we 
study the combined reconstruction from multiple chopped 
data with different chopping amplitudes and the effect on 
the quality of the reconstruction. Finally, in Sect. 0] we 
present some numerical experiments. 



2 



F. Lenzen et al.: Robust Reconstruction from Chopped and Nodded Images 



2. Problem description and basic reconstruction 
methods 

2.1. Continuous problem 

The prob lem description basically follows iBertero et alJ 
l)200 3a*b') . but differs in the way that we take into 
account two-dimension al chopping amplitudes, while in 
iBertero et aP l|2003ald ) the data are preprocessed by ap- 
propriate rotations so that the chopping amplitude is 
along the principal axis of the data. 

We denote by u : R-^ ^ K the brightness intensity dis- 
tribution in the sky, which we assume to be non-negative. 

Let = [0, Ix] X [0, ly] C be a section of the sky 
under investigation and let h — {hx, hy) 7^ 0, hx> 0, hy > 
0, be the chopping amplitude. We define the operator 

u Ih{u){x,y) := 2u{x,y) - u{x + hx,y + hy) (1) 
-u{x - hx,y - hy), (2) 

where Y :— L^{il) and X := L'^{flh) are the spaces of 
square integrable functions on Q and 

flh ■■= i^U {{x,y)\ {x + hx,y + hy) e OT 

{x - hx,y - hy) e ri}, 

respectively. 

The problem of reconstruction from chopped data 
/(x, y) : — > M can be written as an operator equation 

hiu) = /. (3) 

A solution of (O is also a minimizcr of the functional 

u ^ Whiu) - , 

and therefore satisfies the corresponding optimality con- 
dition 

llihu-f) = 0. (4) 
II is the adjoint operator to Ih, which satisfies 

Ih{u)v = / ul^iv) for aU u G L^{nh),v G L^{n) . 

Since is injective, for any injective operator J : X Y, 
©, Q and 

J{hu - /) = (5) 

have the same solutions. Basic numerical methods for solv- 
ing (01 are based on fixed point iterations: 







(6) 



where t is a relaxation parameter. For J = the fixed 
point iteration is commonly referred to as Landweber it- 
eration 



Tl*ih u ^-f) u':=0. (7) 

In IBertero et al.l l)2003a|) . © has been implemented 
with an opera tor J that extends u p eriodically to ^Ih . 
Here, as in IBertero et al.l l|2003bj) . we use the operator 

r, \f \ f u(x) if a; S 
-^("^^'^ ■= I else ' 

which extends a function u E L^{il) to zero in fi/i \ il. 



2.2. Discretization 

The domain CI is discretized by a quadratic grid with nodes 
{xi,yj) e fl,i ^ 1 . . .N,j = 1 . . .M and cell length 1. The 
nodes coincide with the sampling points of u and /, i.e., 
u := (mjj) = {u{xi,yi)) and / := (/ij) are the sampling 
data. 

With u a bilinear interpolating function is associated. 
Assuming w = on fi/i \ the resulting discretized system 
of ® is 



AhU = f . 



(8) 



where Ah is an NM x NM matrix. Details on the struc- 
ture of m atrix Ah and properti es of its eigenvalues can be 
found in IBertero et alJ l|2003bj) . Note that matrix Ah is 
symmetric and positive definite and thus (jSJ has a unique 
solution. 

Since the null-space of Ih is not trivial, the choice of 
the extension of m into ^h\^ enforces a particular solution 
of Ih to be calculated. 

2.3. Review on reconstruction methods 

Three different methods have been proposed in the lit- 
erature for reconstruction from chopped and nodded 
data: Fourier-base d reconstruction method (cf. iBeckersI 
[Bgr^groe^^^ 2003a') . iterative reconstructions (cf. 
Bertero et al.ll2nn3al) . and a wavelet based approach (cf. 
Chan et alJl2003f) . These methods are reviewed below: 



1. For applying the F ourier-based reconstruction (cf. 
IBertero et al.l (12003a') 'I it is assumed that h is in verti- 
cal direction and the solution of ^ is periodic across 
dCl in chopping direction. In this case the image can 
be reconstructed for each column separately, the dis- 
cretized linear system of Ih becomes AhU = f with 
a circulant matrix Ah and Fourier techniques can be 
used for the efficient numerical solution: applying the 
discrete Fourier transform to the linear system it be- 
comes 

HAh):Fu = Tf 

where ^F{Ah) = D := diag{di, . . . dn) is diagonal, and 
thus can be solved efficiently provided the diagonal 
matrix has full rank. 

2. The c onstrained Landweber iteration ( Bertero et alJ 
2003a) (referred to as method (A)) reads as follows: 



(0) 







u 

tor n — 1 ... N 

= P+ 



An-l) 



tAI {f-Ahu(^' 



-1) 



)] 



where P+ is the projection operator onto the set of 
non-negative vectors (each negative entry is set to 
zero), and r is a positive relaxation parameter. 
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The projected Lavrentiev iteration (method (B)) is de- 
fined by 

7.(0) = 

for n = 1 . . . iV : 

Both iterative methods work reasonably efficient if the 
chopped and nodded data are only distorted by a small 
amount of noise. In this case, the results are qualita- 
tively comparable to those obtained with the Fourier 
methods. However, these methods suffer from robust- 
ness with respect to high noise distortions. 
3. Finally we review a wavelet approach proposed by 
IChan et~al\ l)2003|) . Here chopping is interpreted as 
high-pass filtering. Supplementing this filter by a low- 
pass and a subsequent high-pass filtering a tight frame 
wavelet system for a multi-resolution analysis is ob- 
tained. The Landweber method is then combined in 
a multi resolution framework and wavelet threshold- 
ing is applied (cf. Donoho & Johnstone 1994). Wavelet 
thresholding is used for denoising in each iteration 
step. Moreover, a post-processing step for removing 
artifacts in feature-less areas is proposed. 

3. Robust reconstructions 

3.1. Modifications of the iterative metliods 

Simulations with noisy data show that the methods de- 
scribed above have a tendency to introduce artificial struc- 
tures from the noise (cf. Sect. I4.2jl . 

We propose to combine iterative reconstruction meth- 
ods with a median filtering tech nique in each i teration 
step. The median filtering (cf. e.g. IPestmaiJIlQQSj) ) is the 
method of c hoice, since it removes artifical structures (cf. 
ISoil]pll2nn,'^ appearing in each iteration step. Moreover, 
the median can be implemented very efficiently. 

The median filter is defined as follows: 

For an odd number of values {vi . . . V2n+i} in ascend- 
ing order, the median is w„. In median filtering each value 

jp is replaced by the median value of surrounding val- 
ues Uij with indices i,j in a neighborhood of io,jo. We 
take as neighborhood («, j) satisfying \i — io\ < s and 
li^iol s with s = 1, 2, 3. The median fi lter c an be imple- 
mente d very efficiently (cf. ISoilld l)2003(l and ISonka et al.1 

We investigate the following variant of method (A) 
defined by 

= 

for n^l...N : 

1 . = P+ (u("-i^ +TAl{f- 

2. Apply the median filter to u'"-* . 

and the variant of method (B) is defined accordingly. 

This strategy of combination of an iterative method 
with additional filtering after each iteration is analogous to 
the strategy of the wavelet based approach in lChan et alJ 



l|2003^ . where wavelet thresholding is used for additional 
filtering after each iteration of the Landweber method. 

3.2. Reconstruction based on the conjugate gradient 
method 

Since Ah is symmetric and positive-definite, it can be 
solved with the conju gate gradient (eg) method (cf. 
iHanke-Bourgeoij 120021) . which has faster convergence 
properties than the Landweber and Lavrentiev methods. 
However, the fast convergences also makes the method 
more sensitive to noise, which we overcome again by ap- 
plying after one iteration of the cg-method an additional 
median filtering step. The modified cg-method reads as 
follows: 
w(0) = 

for n^l...N : 

1 . Apply one step of the eg -method with 
initial vector u^"^^' 

2. Denoting the result by u^") , 

we replace u*-") by P+u*-") to meet the 
constraint of non-negativity. 

3. We apply the median filter to {t*^") 
to achieve iterate u*-"' . 

The constraint of non-negativity is necessary to avoid 
artifacts in the reconstruction. 

3.3. l\/lultiple chopped data sets 

For an improvement of the quality of the reconstruction 
lBecker£ll)l994l) :l Bertero et al.l ()2000|l proposed to use mul- 
tiple chopped and nodded data sets. The reconstruction 
is performed on each data set independently and the re- 
sults are combined in a post-processing step calculating 
the p ointwise mean resp. the median. 

In iBertero et al.l (jipOO) this strategy is used mainly 
to avoid two kinds of artifacts, first, ghosts from bright 
sources, and second regions, where faint structures are su- 
perimposed by negative counterparts of bright sources and 
consequently the reconstructed image becomes zero artif- 
ically. 

In the following we take advantage of multiple data 
sets for a robust reconstruction in the presence of noise. 

Let h'^jk = 1 . . . K denote a set of chopping amplitudes 
and denote by f'' the corresponding sampled data sets. 

A reconstruction w is a solution of the system 

Ah.u^ f k^l,...,K . (9) 

We solve this system with a blocked Landweber- 
Kacmarcz method ( for m ore details on such methods cf. 
iKaltenbacher et all 1)20051) ): 

Set = 

for n=l...N: (iterationindex) 

1 . For each fc = 1, . . . , A' perform 
two eg iterations starting from 
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(u("-i)) . 

(n) 

The solution is denoted by uj^ ■ 
2. Calculate the median u^"-* from 
of u^i^\ . . . , u^^"^ for noise removal. 

Note that the median is calculated separately in each 
sampling point for the stack of images: That is in step 2 
of the above algorithm we set 

:= median {(ui)ij, . . . , {uK)i,j} 

for each i, j. 

We also tested an alternative combination, namely to 
first detect outliers in {(ui)i.j, . . . , {uK)i,j} via the stan- 
dard deviation and then calculating the mean ignoring 
these outliers, which gives fairly similar results as when 
using the median. 

4. Numerical Experiments 

4.1. Simulation of the Noise Process 

Realistic noisy test data, which are used for our numeri- 
cal experiments, require detailed knowledge on the noise 
process. There are two different sources of noise, which 
we have to be taking into account, background emission 
and thermal (dark) noise which is aimed to be reduced 
with chopping and nodding, and noise occurring during 
the data recording process, which is Poisson noise (due to 
photon counting errors of a CCD detector) and Gaussian 
read out noise. 

We assume that the noise affects each of the four 
recorded signals independently. Let ri{u{x,y)) denote the 
random variable, then the recorded signal reads as follows 

J{x, y) := 2 {u{x, y)) + ry(u(a:, y)) + i]{u{x, y)) 

-u{x + K,y + hy) - ri{u{x + h^.y + hy)) (10) 
-u{x -h^,y~ hy) - ri{u{x - lix,y- hy)) . 

We write f{x, y) := f{x, y) + f]{x, y) resp. f,j := fij +Tiij, 
where fj{x,y) 

fj{x, y) = Tjiuix, y)) + ?7(u(a;, y)) 

-ri{u{x + hh,y + hy)) - ri{u{x - h^^y - hy)). 

According to our considerations, the random variable 
77 (u) is the sum of a Gaussian white random rji = 
iV(0, (Ti), (Ji > (the background emission and read out 
noise) and a Poisson noise ri2{u): ri{u) := 771 (u) + ri2{u). 

We use test images with intensities scaled to the range 
of [0,1]. These intensities are related to the number of 
photons counted by the CCD detector. An intensity u cor- 
responds to a number of N ■ u photon counts with some 
unknown factor A'^ > 0. The Poisson distribution of the 
photon counts can be approximated by a Gaussian distri- 
bution with standard deviation ^/N ■ u. 

Scaling this distribution to the range [0, 1] the inten- 
sities are Gaussian distributed with standard deviation 



y^u/N. Therefore, we can approximate r]2 with a Gaussian 
distribution with standard deviation ^Ja^u and u^ — j^- 

4.2. Behavior of the 'classical' methods 




Fig. 1. Left: Test image. The white frame marks the 
domain il, where the chopped data are recorded; right: 
Corresponding chopped data with chopping amplitude 
h = (5, 3) including noise with variances cri = 0.05 and 
0-2 = 0.0001. 



In the following we demonstrate that artificial 
structures appear in the reconstruction with 'classical' 
methods from the noise. 

Hereby we concentrate on the case of chopping 
amplitudes being small with respect to the image size. 

In the case of a horizontal chopping amplitude being an 
integer multiple of the pixel size the condition number of 
the matrix is of the order of the ratio of n umber of sam- 
pling points and the chopping amplitude fcf. lBertero et all 
2003b). Thus the problem becomes ill-conditioned for 
small chopping amplitudes and low-frequency eigenvectors 
are affected by noise i n the data resulting p roducing the 
artificial structures fcf. iBertero et al]l2003b|) . 

Fig. n shows an artificial test image (left) and corre- 
sponding chopped and nodded data (right) with chopping 
amphtude h — (5, 3) which is distorted 

by two Gaussian noise processes with parameters ai = 

0.05 and 0-2 = 0.0001. The signal-to-noise-ratio is 25.4752 
1 

Fig. 121 top left and top right, resp., shows the recon- 
struction from these test data applying method (B) with 
T — 0.1 for h = (5,3) after 10 and 100 iterations, re- 
spectively. We applied the same method to test data with 
a chopping amplitude h = (5.5, 3.3) not matching the 
grid spacing. Reconstructions after 10 resp. 100 steps are 
shown in Fig.|21 bottom left resp. bottom right. 

^ Let M : n ^ R be a signal with minimum Umin and maxi- 
mum Umax, S : f2 ^ K a distorted signal and a the standard 
deviation of u — u, then the signal-to-noise-ratio is defined by 
SNR := 201ogio( ""°°'r"''" )- 
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The computation times of method B were 0.05 sec- 
onds for h = (5, 3) and 10 steps, 0.2 seconds for h — 
(5.5,3.3) and 10 steps, 0.5 seconds for h = (5,3) and 100 
steps and 1.9 seconds for h — (5.5,3.3) and 100 steps. 

Since the chopped and nodded data / provide informa- 
tion at the objects' edges, their reconstruction is satisfac- 
tory at an early stage of the iteration, whereas the objects' 
interior is reconstructed at a later stage. Therefore, larger 
objects require a larger number of iterations to be recon- 
structed satisfactory. 

When the discretization points are aligned with the 
chopped discretization points the reconstruction is more 
efficient. If the points are not aligned, then the numer- 
ical reconstruction is smoother; moreover, the iterative 
algorithms are slower convergent, i.e. a larger number 
of iterations is needed for reconstructing the interiors of 
the objects, which in turn leads to an amplification of 
noise and appearance of artificial structures. Method (A), 
the cg-method, and the Fourier method shows similar 
convergence properties. Applying the cg-based method to 
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Fig. 2. Top left: Reconstruction from the noisy data given 
in Fig. n right, with method (B) after 10 iterations and 
r = 0.1. Top right: Reconstruction with method (B) after 
100 iterations and r = 0.1. Bottom left: Reconstruction 
with h = (5.5,3.3) using method (B) after 10 iterations 
and T ~ 0.1. Bottom right: Reconstruction from the same 
data using Method (B) after 100 iterations and r = 0.1. 



noisy data shows that the residual \\AfiU — f\\ decreases 
in the beginning but after some iterations is starts 
oscillating or is even increasing. This suggests to stop the 
iteration of the cg-solver when the residual reaches its 



first local minimum. 

Since Eq. (PJ is linear we have Ih{uf -t- Wj;) = / + ?7, 
where w,, are the data reconstructed from noise and u/ is 
the reconstruction from noise free data. 

For illustration of the noise amplification behavior of 
this equation, we assume that the data / = are disturbed 
at one grid point: 



1 if i = io and j 
"^'^^ = i else 



Jo 



for a fixed io £ {I . . . iV}, jo e {1 • • • M}. 

First, let us consider the Fourier-based reconstruction 
from this particular 77 for horizontal chopping amplitudes. 
In order to handle small eigenvalues of matrix D, we use 
a slightly modified method replacing matrix D by 

Dl := diag{d\ ... 4) 

where £ > is a small parameter and 



d] 



l/d» 



if \di\ > £ 
else 



1 . 



Fig. U] shows a horizontal cut through point {io,jo) 
of the function reconstructed from 77 for h = (5, 0) and 
£ = 0.1 resp. £ = 1. 




Fig. 3. Horizontal cut of the image reconstructed with 
the Fourier method with Dq i (solid line) and Di (dashed 
line). The data are fi,j := Si„idj„j for afixed io S 
{1 . . . N},jo E {1 . . . M} and the chopping amplitude is 
h = (5,0). 



Similar artificial structures can be observed in the re- 
construction from iterative methods. Fig. 0] shows the 
reconstruction from data / for h = (5, 0) (top) and 
h = (5.5,0) (bottom). 

In the later the unknowns of u are far more 

coupled the reconstruction is smooth and thus more com- 
plicated to handle numerically. 



^ Computed on an AMD 64 FX 3500-I-, computations times 
have been averaged over several runs 
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Fig. 4. Reconstruction from data = {1 if i ~ 
io and j — jo,0 else }, io and jo fixed, using the cg-based 
method for h ~ (5,0) (top) and h = (5.5,0) (bottom). 



4.3. Results of the modified methods 

In this section we present numerical results of the modified 
methods applied on noisy data. 

Firstly we used the Fourier method with inverse to 
reconstruct our first test image. ^ 

But numerical experiments show that a small e is 
needed for a good reconstruction for noise free data. Since 
the eigenvalues of the inverse tend to infinity for e ^ 0, 
artificial structures due to noise are significant for small 
e. For noisy data, a compromise between good reconstruc- 
tion and little noise amplification is not possible. The 
results are not satisfactory, so they are not presented here. 




Fig. 5. Left: Reconstruction from the data in Fig.^ right, 
with chopping amplitude h = (5, 3) using method (B) with 
T = 0.1 and 100 iterations combined with median filtering 
(filter size s = 1). Right: Reconstruction after 30 iterations 
of the cg-based method combined with median filtering 
(filter size s = 3 resp. s = 1 in the last iteration). 



Secondly, we consider reconstructions from the cg- 
based method and method (B). Results from method (A), 
which in general look more blurry, are not presented here. 

Fig. El left, shows the reconstruction from noisy data 
with chopping amplitude h — (5, 3) with method (B) com- 
bined with median filtering (filter size s = 1). Fig.|Slright, 
shows the reconstruction using the cg-method with me- 
dian filter (s = 3 resp. s = 1 in the last step). 



The computation times were 3.3 seconds for method 
(B) and 4.4 seconds for the cg-method. 

In both reconstructions artificial structures arising 
from noise are removed by the median filtering (cf. Fig. 12). 




Fig. 6. Left: The second test image. Right: the corre- 
sponding noisy chopped and nodded image for h — (10, 3), 
CTi = 0.002 and era = 0.00001. 




^ For implemen tation the FFTW-C-library 

l|Frigo fc ,Tohnsonlll99Sl. resp. www.fttw.org) was used. 



Fig. 7. Top left: Reconstruction with the cg-method af- 
ter 40 iterations (without median filtering). The chopping 
amplitude is h = (7, 14). Top right: 40 iterations with the 
cg-method combined with median filtering with filter size 
s = 3 (s = 1 in the last iteration). Note that for com- 
parison with the method of multiple chopped data sets 
with amplitudes h = (10,3), (0,7), (10,10), (7,14) and 
(15, 0) we depicted h = (7, 14) showing best results by in- 
dependent reconstruction. Bottom left: Combined image 
using the multiple chopped data calculating the pointwise 
median after independent reconstruction. Bottom right: 
Reconstruction using multiple data sets and combining 
the results after each iteration using the pointwise me- 
dian. 
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Further experiments have been performed with the test 
image shown in Fig. which contains an object with a 
wide halo extending across the boundary of the chosen 
domain. 

Fig. El (left) shows the signal u on flh, the part in which 
the chopped data are provided is marked by the white rect- 
angle. We produced a multiple chopped data set with 5 dif- 
ferent chopping amplitudes (10,3), (0,7), (10,10), (7,14) 
and (15,0). The chopped data for h = (10,3) on il is 
shown in Fig. |H1 (right). 

We added noise of variance <Ji — 0.001 and (72 = 
0.00001. Note that for this test image the chopped data set 
contains very weak structures. The signal-to-noise-ratio is 
40.0494. We compare different strategies for reconstruc- 
tion from noisy data. First of all we calculated the re- 
sults of the cg-method with 40 iterations for each single 
data set, applying the median filter after each iteration. 
For a fair comparison we depict the data set showing the 
best results, which is that for h = (7,14). Fig. [7| (top 
left) shows the result of the cg-based method without 
median filtering. Undesirable noise enhancement occurs. 
Fig. [71 (top right) shows the reconstruction with the cg- 
method where in each iteration step a median filtering is 
performed (s = 3 resp. s = 1 in the last step). 

Fig. [3 (bottom left) shows the combined result of inde- 
pendent reconstructions from the multiple data set, apply- 
ing 40 iterations of the cg-method without median filtering 
on each set and finally calculating the median in each data 
point. By this post-processing the artificial structures are 
weakened. 

Fig. [71 (bottom right) shows the reconstruction with 
combining the different results after each cycle of cg- 
iterations by calculating the pointwise median. 

The computation times were about 0.9 seconds for the 
cg-method without median filtering and about 5.5 seconds 
with median filtering. The reconstruction with multiple 
chopped data sets took about 4.7 seconds with indepen- 
dent reconstruction and finally averaging and 5 seconds 
with calculating the median after each iteration. 




Fig. 8. Left: Third test image, planetary nebula Menzel 
III, Right: Artificial chopped data for h — (14, 10) includ- 
ing noise ci = 5 10~^ and 1T2 = 10^^. 



A third test image (Fig. [SJ is an observation of the 
planetary nebula Menzel III (the ant nebula). The inten- 
sities stored in floating point precision show a high dif- 
ference of intensities between the nebula and the central 
star. To visualize the data, the intensities were logarith- 
mically scaled to gray values by applying the function 

Six) = log{x/{rn-raax)+OA)-log{a.l) ^ > q . 

^ ' log{lA)—log{Q.l) ' — ' 

the maximal intensity in the image and m is a scaling 
parameter, set to m = 0.01 for the results presented here. 

We simulated the chopping and nodding procedure 
including artificial noise of variance ai = 5 • 10"^ and 
(72 = 10~^. We have chosen this noise variance correspond- 
ing to the weak intensities of the nebula in the range about 
10~^. We used small chopping amplitudes up to 14 pixels 
for testing. Note that when applying large chopping am- 
plitudes the original image can be reconstructed by a few 
iteration steps. Then the amplification of noise is negligi- 
ble. The corresponding chopped data for h = (14, 10) is 
given in Fig. [HI (right). The signal-to-noise ratio is 42.9788. 
Fig.[ni(top left) shows the result of method (B) after 100 



Fig. 9. Top left: Result of method (B) after 100 iterations 
in combination with the median filtering with filter size 
s = 1. Top right: Result of the cg-based method after 20 
iterations in combination with median filtering with size 3 
resp. 1 in the last iteration. Bottom left: Result of a com- 
bined reconstruction with method (B) (r — 0.1, 100 steps) 
using multiple data from chopping amplitudes (14,10), 
(14,0), (14,14), (10,0), and (7,14). After each iteration 
step the reconstructed images are reinitialized with the 
pointwise median. Bottom right: Result of a combined re- 
construction (pointwise median) with the cg-method with 
20 steps using the same multiple data set. 



iterations where at each iteration median filtering with 
s = 1 has been used. Fig. [HI (top right) shows the result 
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of the cg-based method after 20 iterations in combina- 
tion with median filtering. In both cases some blurring ef- 
fect caused by the median filter can be recognized. Finally 
we used data obtained with chopping amplitudes (14, 10), 
(14,0), (14,14), (10,0) and (10,14). Fig. H (bottom left) 
shows the reconstruction using method (B) on each data 
set and combining the results after each iteration step us- 
ing the pointwise median as described in Sect. 13.31 In 
Fig. El (bottom right) we used the cg-method (20 itera- 
tions) where again the pointwise median is calculated after 
each iteration . 

The computation times were about 2.5 seconds using 
method (B) with median filtering, 0.9 seconds using the 
cg-method with median filtering and 2.8 resp. 1.6 seconds 
using method (B) resp. the cg-based method on the five 
different chopped data sets and calculating the median 
after each iteration. 



5. Summary and Conclusions 

In this paper we give an overview of Fourier-based and it- 
erative schemes for reconstruction of images from chopped 
and nodded data. As an alternative to reconstruction 
methods documented in the literature we proposed an 
algorithm based on the conjugate-gradient (cg)-mcthod 
which is faster converging. 

Experiments show that for the 'classical' reconstruc- 
tion methods noise can severely affect the reconstruction 
process. To prevent noise enhancement we propose to com- 
bine iterative methods with median filtering techniques 
(robust methods from statistics). Filtering is performed 
in each iteration of the reconstruction process. 

We produced satisfactory numerical results for three 
different test images with noisy chopping and nodding 
data. Noise artifacts are removed during reconstruction. 

An alternative strategy to enhance the quality of re- 
construction is to use multiple chopping amplitudes. The 
results of the reconstruction are combined by statistical 
methods after each iteration. This also provides a method 
robust against noise. 

Since the method of multiple chopped data sets does 
not show any blurring effects, we propose to use this 
method if the required data are available. In the case of a 
single data set the methods with median filter should be 
applied. 
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